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Abstract. We review the application of perturbative QCD to diffractive processes. 
We introduce the two gluon exchange model to describe diffractive qq and qqg 
production in deep inelastic scattering. We study the triple Regge limit and briefly 
consider multiple gluon exchange. We discuss diffractive vector meson production 
at HERA both at t = and large \t\. We demonstrate the non-factorization of 
diffractive processes at hadron colliders. 



1. Introduction 

Diffractive scattering has become a vast area of study in particle physics and generated 
a wide range of theoretical approaches. Here we will concentrate on those processes 
which have the potential to provide a deeper insight into the dynamics of strong 
interactions, that is into quantum chromodynamics (QCD). Our main focus will be on 
diffractive deep inelastic scattering. Diffraction in pp-scattering or photoproduction 
will be discussed fairly briefly. The motivation for studying diffractive deep inelastic 
scattering is not only driven by the fact that we have excellent data, but also 
that it offers an opportunity to explore the interesting transition from "hard" to 
"soft" physics. Hard physics is associated with the well established parton picture 
and perturbative QCD, and is applicable to processes for which a large scale is 
present. Soft dynamics on the other hand, associated for example with total cross 
sections and elastic scattering of hadrons, is described by non-perturbative aspects of 
QCD. Diffractive deep inelastic scattering offers the opportunity to directly probe the 
semihard transition region and to link the, otherwise distinct, regimes of soft and hard 
physics. The means by which we attack the problem is perturbative QCD carefully 
extrapolated into the semihard regime. It is thus natural to put strong emphasis on 
diffractive processes, which are believed to be purely perturbative, such as diffractive 
vector meson production in deep inelastic scattering or in photoproduction when the 
momentum transfer is large. 

The classical definition of diffraction in hadron-hadron or (virtual) photon-hadron 
scattering is the quasi elastic scattering of one hadron combined with the dissociation 
of the second hadron or photor||. For example, at the electron-proton collider HERA, 
diffractive deep inelastic scattering is usually written 

7* +p ^ X + Y (1) 



I We consider here only the diffractive dissociation of one of the incoming particles. 
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where Y denotes the elasticahy scattered proton (or one of its low mass excitations) 
and X denotes ah the hadrons coming from the dissociating virtual photon. The 
classical method of detection is to tag the elastically scattered hadron. In the present 
high energy collider experiments, however, the elastically scattered protons usually 
disappear into the beampipe requiring special detectors far down, and adjacent to, the 
beampipe for their detection. Very often these detectors cannot deliver high enough 
statistics or good enough energy resolution so that an alternative method, the rapidity 
gap method, is employed. The rapidity gap method is based on the simple fact that 
the elastically scattered proton (or Y state) travels in roughly the same direction 
as the original hadron leaving a large gap between the rapidity of this particle and 
the other outgoing hadrons forming the inclusive state X. Rapidity is a measure of 
the component of a particle's velocity along the incoming proton direction. That is 
diffraction at the HERA electron-proton collider is characterized by the quasi-elastic 
scattering of the proton and the dissociation of the virtual photon. Such a diffractive 
event is in contrast to the usual deep inelastic scattering event 

1*+P^X (2) 

in which the virtual photon breaks up the proton and the outgoing hadrons populate 
the full interval of rapidity. Nevertheless about 10% of all events in deep inelastic 
scattering are of diffractive nature and high statistics data samples have been collected 
at HERA §. 

Before this diffraction was observed in hadron collisions and is still being measured 
at the Tevatron proton-antiproton collider at the present time. The relative rate of 
diffraction compared to the total rate is on a similar level to that of deep inelastic 
scattering The problem, however, is that the total cross section is huge and only a 
tiny fraction of events are of special interest such as heavy particles {W and Z bosons 
or the top quark) and jets of high transverse momentum. Since the hadronic activity 
is much higher in hadron-hadron scattering compared to deep inelastic scattering the 
detection of rapidity gaps is a tedious task. With a trigger set on W bosons or high-p^ 
jets the fraction of events with a rapidity gap is only on the level of 1% Q. It requires 
a special effort to extract diffractive processes at the Tevatron whereas at HERA they 
almost come for free. 

There is a group of processes which do not precisely belong to diffraction in the 
way we have defined it. These are processes with large momentum transfer across 
the rapidity gap. Examples are the production of opposite di-jets with a rapidity gap 
in between them at the Tevatron, or diffractive vector meson production at HERA. 
Due to the large momentum transfer the proton (or antiproton) always breaks up. 
Because of the large rapidity gap, however, they are closely related to the traditional 
diffractive processes. 

The early attempts to describe diffractive scattering were based on Regge 
phenomenology where the Pomeron is thought to be the leading Regge pole with a well 
defined and unique, process-independent Regge trajectory. The Pomeron trajectory, 
together with secondary trajectories, allows one to fit all kinds of hadronic data (total 
and elastic cross sections) [Q. These universal trajectories anit) form the backbone 
of Regge phenomenology. According to early analyses it seemed that diffraction can 
be described in terms of triple Regge diagrams (see Section 3) without changing the 
trajectories or the couplings to hadrons. The only free parameters were the couplings 
of the Reggeons amongst themselves. The more recent results |^ from the Tevatron, 
however, indicate that this simple picture does not hold at very high energies. The 
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(b) 




Figure 1. Diagrams for DIS, 7*p X, and inclusive diffractive DIS. j*p Xp. 
In the latter process there is a rapidity gap between the outgoing hadrons X and the 
slightly deflected proton. In the diagram the gap has been associated with Pomeron 
exchange carrying a fraction xp of the proton's momentum. 

extrapolation of the cross section from the SppS CERN collider energy to the Tevatron 
energy clearly overshoots the data. The explanation for the discrepancy is most likely 
given by unitarity corrections. 

After introducing the structure functions that can be measured in diffractive 
deep inelastic scattering, wc outline in Section 3 the Rcgge description of diffractive 
processes based on the Pomeron, or the "soft" Pomeron as it is now frequently called. 
The perturbative QCD description of diffractive processes, which is the subject of 
this review, is based on two gluon exchange. This will form the framework of our 
discussions from Section 4 onwards. Two gluon exchange is the simplest form of 
vacuum quantum exchange and, for historical reasons, people often speak of the 
"hard" Pomeron. However the connection between "soft" (non-perturbative) and 
"hard" (perturbative) diffraction is far from clear. The different energy behaviours 
are very evident, as will be seen, for example, when we discuss diffractive vector meson 
photo- and electro-production (see, in particular. Fig. 9). 



2. Diffractive structure functions 

First let us recall (inclusive) deep inelastic electron- proton scattering, ep eX, where 
X represents all the fragments of the proton which has been broken up by the high 
energy electron. The basic subprocess 7*p — > X is shown in Fig. 1(a). It can be 
expressed in terms of two functions F2 and Fl which characterize the structure of 
the proton. These proton structure functions depend on two (invariant) variables, the 
"virtuality" of the photon = —q^ and the Bjorken a;- variable 

x=^- (3) 

where p and q are the four-momenta of the proton and virtual photon respectively, 
see Fig. 1(a). W is the total 7*p centre-of-mass energy. If we were to view the proton 
as made up of massless, point-like quark constituents, then it is easy to show that x 
is the fraction of the proton's momentum carried by the quark struck by the virtual 
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photon. In this simple quark model -F^ = and 

F2 = Ft = (4) 

is independent of Q^. The sum is over the flavours of quarks, with electric charge 
(in units of e) and distributions q{x). Ft,l are the proton structure functions for DIS 
by transversely, longitudinally polarised photons. 

The general form of the DIS cross section, up to target mass corrections, is 

where a is the electromagnetic coupling. The third variable y is needed to fully 
characterize the DIS process, ep — s- eX, namely y = Q^/xs where -^/s is the total 
centre-of-mass energy of the electron-proton collision. 

Now we turn to inclusive diffractive DIS, ep eXp, for which the subprocess is 
7*p Xp, where a slightly deflected proton and the cluster X of outgoing hadrons are 
well separated in rapidity. In terms of Regge terminology the rapidity gap is associated 
with Pomeron (or vacuum quantum number) exchange, shown by the zig-zag line in 
Fig. 1(b). The 'y*p — + Xp subprocess is described by diffractive structure functions 
F^j^, which now depend on four variables: /3 = x/xp, Q^, xp and t. The variable xp 
is the fraction of the proton's momentum that is carried away by the Pomeron. The 
variable /3 in diffractive DIS plays an analogous role to x in DIS 

o _ Q Q fc\ 

' ^ ^ ~ 2xpp.q " + ^ ' 

where M is the invariant mass of the diffractively produced X cluster in Fig. 1(b). 
In analogy to (||), the general form of the diffractive DIS cross section is 

d'^a^{ep-> eXp) Stto^ r \2t TpD{i) m r,^ 

dpdQUxpdt ^^l[^ + (^-^)J^2 iP.Q.^r^t) 

-y^F^'^'^\p,Q^xp,t)]. (7) 

The superscript (4) is to indicate that the structure functions depend on four 
independent variables. The Fl contribution is often neglected due to the smallness 
of y^ in most of the measurements. The variable t — [p — p'Y is the square of the 
4-momentum carried by the Pomeron. Most diffractive DIS events occur for small 
values of t. Usually t is integrated over and measurements of the structure function 



F^''>{P,Q\xp) = / dtF^'''>{p,Q\xp,t) (8) 

are made. 



3. Regge approach to diffraction 

The early attempts to describe diffractive scattering were based on Regge theory (see, 
for example, Q), in which the basic idea is that sequences of hadrons of mass rrii 
and spin ji he on Regge trajectories a{t) such that a{mf) = ji. Prior to QCD, 
strong interactions were thought to be due to the exchange of complete trajectories 
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of particles. Indeed the Regge model is able to successfully describe all kinds of 
"soft" high energy hadronic scattering data: differential, elastic and total cross section 
measurements. In this model the high energy behaviour of a hadron scattering 
amplitude at small angles has the form 

~ ^ (9) 

R 

where for simplicity we have omitted the signature factor, s is the square of the centre- 
of-mass energy and ~t is the square of the four-momentum transfer. The observed 
hadrons were found to lie on trajectories aji{t) which are approximately linear in t 
and parallel to each other. That is hadrons of increasing spin and mass, but with the 
other quantum numbers the same, lie on a single trajectory afj{t). The leading such 
trajectories are the p,a2,oj and / trajectories which are all approximately degenerate 
with 

anit) ~ 0.5 + O.m. (10) 

For example only the p trajectory has the appropriate quantum numbers to be 
exchanged in the process Tr~p — > n'^n. The energy or s dependence of the differential 
cross section da/dt therefore determines ap{t) for t < 0, see (H). For small \t\ the 
trajectory ap{t) is found to be linear in t and, when extrapolated to positive t, to pass 
through the p{l^) and p{'S^) ■ ■ ■ states, i.e. ap(vn?p) = 1, 3, ... at the appropriate mass 
values. 

So far so good, but the total cross sections are observed to increase slowly with 
energy at high energies, so we need a higher lying trajectory. This is best seen from 
the optical theorem which expresses a total cross section (say, for AB ^ X) va terms 
of the imaginary part of the forward {AB AB) elastic scattering amplitude 

a{AB^X) = ilmA(s,0) = V /3j? s"«(°)-^ (11) 

R 

To account for the asymptotic energy dependence of the total cross sections a 
Pomeron^ (or vacuum quantum number exchange) trajectory is invoked with intercept 
aip(O) ~ 1.08. Indeed the total, elastic and differential hadronic cross section data 
are found to be well described (for small \t\) by taking a universal pole form for the 
Pomeron 

ap{t) ~ 1.08 + 0.25i, (12) 

together with the other sub-leading trajectories as in ( |l0|) Q. The Pomeron should be 
regarded as an effective trajectory, since the 5" °* power behaviour of the total cross 
sections will ultimately violate the Froissart bound. The link between this successful 
Regge description of "soft" processes and the underlying fundamental theory of QCD 
is not, as yet, known in detail. Most probably Pomeron exchange originates mainly 
from the exchange of a two-gluon bound state, while the meson trajectories (p, 02, w, /) 
correspond to qq bound states. The Regge Pomeron discussed above is now often called 
the "soft" Pomeron. 

To apply this approach to inclusive DIS and, more especially to its diffractive 
component, we again make use of the optical theorem, together with its generalisation 

t Originally the total cross sections were thought to asymptote to a constant at high energies and 
so a Pomeron with ap (0) = 1 was introduced tel . 
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(a) DIS (optical theorem) Regge limits 

o(fp^X) ~ F2(x,Q') 





Figure 2. The DIS and diffractive DIS cross sections expressed symbolically 
as structure functions via the optical theorem and its generalisation. The Regge 
limits are shown in the rightmost diagrams (where the photon couples via a quark 
configuration). A summation over allowed Regge exchanges is implied in (b). 

by Mueller Q, which are symbolically shown by the equalities in Fig. 2. The optical 
theorems express the total cross sections in terms of the imaginary parts of the 2-body 
(or 3-body) forward elastic scattering amplitudes, or to be precise the discontinuities 
of the amplitudes across the cuts along the (or IvP) axes, which are indicated by 
the dotted lines in Fig. 2(a) (or (b)). The last diagrams show the various Regge limits 
for the structure functions, where the coupling to the photon is via a quark line. For 
DIS this gives 

i i 

for small see (|). In the naive parton model the valence and sea quark contributions 
to F2 are associated with meson and Pomeron exchange respectively, and so using (^ 
we have 



xqv ~ a;i-"«(°) - x' 
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(14) 

for small x. 

For difFractive DIS, j*p Xp, we apply Mueller's generalisation of the optical 
theorem [Q. For this difFractive case, that is when s/M^ is large, the theorem is shown 
pictorially by the first two diagrams of Fig. 2(b). That is the cross section is given 
by the discontinuity across the cut of the (three-body) 7*pp elastic amplitude, 
where a sum over the exchange Reggeons is implied. The Regge prediction depends on 
whether is large or small. For small the quark box gives the main contribution 
to photon-Pomeron scattering. In Ref. jsj a charge conjugation C = 1 vector current 
coupling of the Pomeron to quarks was introduced with a free coupling constant. The 
resulting contribution to diffraction due to the quark box diagram is found to be 

F,^ ^ (3(1 (15) 

This model shows already many of the properties which are also present in the two 
gluon exchange approach |^ ^ It is important to note the leading twist nature of 
diffractive deep inelastic scattering, that is F.P is independent of Q^. One also finds 
the strong alignment of the quark-antiquark pair along the photon-Pomeron axis which 
is a consequence of the low transverse momentum fcj of the quarks relative to this 
axis. The particular property of alignment has already been pointed out by Bjorken 
Ip^ . It was argued in Ref. Q that the quark-Pomeron coupling for inclusive DIS and 
diffractive DIS should be the same and by fitting inclusive data one can extract the 
coupling constant and make predictions for diffraction. Indeed some 10% of HERA 
DIS events were predicted to be diffractive. 

For large M^, on the other hand, we have the double Regge limit [s/A'P — > oo 
and oo) and the diffractive structure function is described by a sum of triple 

Regge diagrams 

-H^^^^ii^) {Mr^''\ (16) 

The leading behaviour, which is given by the triple Pomeron contribution, is 

F^ - (■^j2^ar(0)-2ajp(t) _ I / . (17) 



xqs 



'(0) 



-0.08 



4. Two-gluon exchange model of diffraction 

We now use perturbative QCD to describe the diffractive DIS process 7*p Xp. In 
QCD the "Pomeron" (or vacuum quantum number exchange) is, in its simplest form, 
represented by two gluons ^] . The minimum number of gluons to form a colourless 
state is of course two. It is not excluded that more than two gluons are exchanged and 
it is important that whenever we talk about two-gluon exchange to remember there is 
the possibility to extend the formalism to multigluon exchange. One might object that 
the whole process is soft and perturbation theory not applicable. Saturation effects for 
high parton densities, however, screen soft contributions, so that a fairly large fraction 
of the cross section is hard and therefore eligible for a perturbative treatment |l^ . 

t The problem that this model breaks electromagnetic gauge invariance, and a possible remedy, has 
been discussed in Ref. P|. 
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The two basic pQCD diagrams are shown in Fig. 3, in which the photon dissociates 
into either a qq pair or qqg. In the first diagram the q and q carry fractions a and 
(1 — a) of the momentum of the photon, and have transverse momenta ±fcf. In the 
second diagram these variables apply to the gluon and the (qq) system. At first sight 
it might appear that qqg would be a small correction to qq production on account of 
an extra as factor. However qqg production dominates at large M^. It incorporates 
an extra t-channel spin 1 gluon in contrast to the (lower lying) t-channel spin i quark 
exchange of qq production. In fact qqg is described by the triple Pomeron diagram 
and gives the leading contribution at large (or small /3) , whereas qq production is 
leading at small M^, see Fig. 2(b). 




Figure 3. Diffractive qq and qqg electroproduction mediated by the exchange of two 
gluons with transverse momenta Ht ■ We view the photon as dissociating into either 
a quark dipole or an effective gluon dipole, made up of a gluon and a compact qq 
pair. 

The intuitive picture is as follows. In the target (proton) restframe the photon 
dissociates into a 55-pair far upstream the target. The (75-pair may radiate a gluon, 
and the whole parton configuration scatters quasi-elastically off the proton via a two- 
gluon exchange. The timescale on which the fluctuation occurs is proportional to 
l/(a; mp) where rup is the proton rest mass. At very small x the fluctuation is long 
lived whereas the scattering is a sudden short impact of the qq-pair or the ggg-final 
state on the target. The impact changes the virtual into a real state but it does not 
change the position in impact parameter space which can be viewed as being frozen 
during the scattering. 

It is informative to show why the fluctuation and interaction timescales are so 
different. For this it is convenient to use light-cone perturbation theory (see, for 
example, [p^) and to express the particle four momenta in the form 

fc^ = {k+,k-,kt) (18) 

where k± — fco ± fca. In this approach all the particles are on mass-shell 

= - k^ = (19) 

and fc+ and kt are conserved at each vertex. In the proton rest frame we have 

= (mp,r7ip,0), = (17+, -Q'^/q+, 0), 

(20) 

k^ = (fc+, r74/fc+, kt), £^ = e^je+, It), 
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with m|i = 
principle the 7 

1 

^ AE ^ 



where m is the mass of the quark. According to the uncertainty 
qq fluctuation time 



g_ - /ci- - k2- 



2q^ 



2q+ 

12 



+ mip/a + m|,/(l — a) 

at high Q^, where /ci and fc2 are the four momenta of the quark and antiquark. The 
momentum fractions, shown in Fig. 3, are a = fci+/g+ and 1 — a = k2+/q+. The 
factor 2 in (21) occurs because the energy ko = + fc_)/2. An estimate of the 
interaction time can be obtained from the typical time for the emission of a gluon of 
momentum £ from the quark ki, say. Then 



- k[_ -t- 



2q+ 



;/a — m^ja' — t1/{a — a') 



2^4 



1 



TTlpXe 



(22) 



since a — a' — £+/q+ and a' = k'i_^_/q^. We can regard xi as the Bjorken variable 
for the gluon-proton interaction. We are concerned with the kinematic region x <^ Xi 
where the leading log(a;f/a;) approximation is appropriate, and so we have ^ n. 

One should note that the above picture is valid in a certain frame in combination 
with a certain physical gauge condition, p ■ A = {p is the proton momentum and A 
the gluon vector potential). In this gauge the parton shower evolves from the photon 
end while Bremsstrahlung from the proton is suppressed. One can turn to the Breit 
frame where the proton is fast moving and with it the Pomeron. The Pomcron in 
this frame is a long lived and very complicated virtual fluctuation. It is then more 
appropriate to use the gauge condition q' ■ A ~ {q' — q + xp is the corresponding 
light-cone vector in the photon direction) under which the parton shower evolves from 
the Pomeron end. The results do not change but the physical interpretation is less 
appealing. 

The virtual fluctuation of the photon into a gg-pair or a qqg-fmal state is described 
by photon wave functions. The QCD-wave function for a qq-pair has first been 
discussed in Refs. 
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Similar results have already been derived within QED. 
In Ref. [ [17[ the fluctuation of a photon into a pair of muons has been studied as a 
model for deep inelastic scattering. The discussion of gluon radiation is a more recent 
development starting with Ref. . 

The wave function for a gg-pair depends on the polarization of the parent photon, 
whether it is transversely (T) or longitudinally (L) polarized. The key assumption 
in the whole approach is the eikonal-type coupling of the t channel gluons with 
momentum £ to fast moving s channel quarks and gluons such that 



{p~£\r\p) - 



(23) 



where p and p ~ £ are the momenta of the fast incoming and outgoing partons. This 
approximation is correct at very high energies {£ <^ p). The derivation of the transverse 
and longitudinal wave functions, based on an explicit spinor representation, can be 
found in Ref. [|9| or Ref. pOl. For a transversely polarized photon of, say, helicity 



A. 



-1 the wave functions V^aa' are 



V2eqakt ■ e{Xj — 1) 
fc2 + a(l-a)Q2 ' 

_ %/2eg(l - a)fcf ■ £(A^ = 1) 
fc2 + a(l - a)Q2 



(24) 
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where e is the polarization vector of the virtual gluon, A, A' are the helicities of 
the quark, antiquark (with ± denoting helicity ±i), and Cq is the electric charge 
of the quark. The general structure of (|2^) is clear. For a massless quark, helicity 
conservation at a quark-photon vertex requires A 7^ A'. The momentum fractions 
a, (1 — a) appear in the numerator because the photon helicity tends to follow the fast 
quark, antiquark. Angular momentum conservation forbids coUinear q, q production. 
Hence the kt factor signifying a P-wave interaction. Finally, the structure of the 
denominator has its origin in ( |2l] ) with rn^ = kf. The longitudinal wave functions 
(with A-y = 0) are found to be 

r r 2e„a(l — a)Q 

The effective wave function for the ggf;-state requires some remarks. Since the 
photon does not couple to gluons directly it cannot decay into a gluon dipole. At 
large Q^, however, the colour structure of the qq effectively combines into that of a 
gluon leading to a gluon dipole. This feature is based on the strong ordering of the 
transverse momenta in the leading \og{Q^) approximation to diffractive DIS, that is, 
it is valid when kt ^ Q and a <C 1, and when the photons are transversely polarized. 
The contribution due to the quark-box can be factorized (as becomes more explicit 
later on) and we are left with the following expression for the effective gluon dipole 
wave function describing the 7* {qq)g transition [pTj |22|: 



1 k1 (£1 -£2) - 2(fct ■ £i)(fct • £2) 



*999 — 1,2 , „^2 ' y-^^) 



where £1 and £2 are the polarization vectors of the (qq) and g forming the gluon dipole. 
Since a <C 1 we have replaced^ (1 — a) by 1 in (|2^). Otherwise the wave function is 
to be interpreted in analogy to eqs. (|2J, |25| ). 

In order to ensure gauge invariance we have to consider all possible couplings 
of the two t channel gluons to the quark- or 'gluon'-dipole. In total there are four 
configurations of which only one is shown in Fig. ^. This leads to the following 
expression 

D-^{kt,£t) = 2^{a,kt) - ^{a,kt+£t) - ^(a,kt-£t) 

a^^^ for^,«fc„ 
* * dkldki (27) 



2*(a,fct) for^t>fc; 



The amplitude for diffractive scattering is basically obtained by folding D'i> with the 
unintegrated gluon distribution T: 

^""^ f^D^{kt,et)as{ft)H^]P,ft) ■ (28) 

The form for is an example of the so-called 'fct-factorization theorem' |Q, here 
applied to the transverse momentum it. The distribution J-', unintegrated over 
contains all the details of the coupling of the two t channel gluons to the proton 

t In principle, there is a symmetry between a and (1 — a), however, by convention, a is used to 
denote the small component. 
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as indicated by the bubble in Fig. This gluon distribution T{x,(^) is a universal 
function applicable to all hard scattering processes involving the proton and, indeed, 
the gluon is the dominant parton distribution at small x. Integrating !F{x,t1) over 11 
gives the conventional gluon distribution xg(x). 

We have now, in principle, all the "QCD-based" ingredients to write down the 
difFractive cross section. Before we do that we comment on the role of the wave 
function for the /3-spectrum. It has been demonstrated in Ref. |Q that the basic 
shape of the /^-spectrum arises from the wave functions. It depends only weakly on 
the details of the unintegrated distribution J^. Looking back at (|2^) we see a 'hard' 
{It <C fct) and a 'soft' {It >• fct)|| limit for . The typical scale for it is given by the 
dynamics inside the proton. The hard limit is represented by the second derivative 
of the wave functions and the soft limit by the wave function itself. The second 
derivatives are 



dk\dkl [/c? + a(l-a)Q2]3 V ' 



^ - 4 e, «(1 - «) Q g-^a-")^; (30) 
dk\dk\ ' ^ [fc2 + a(l -a)Q2]3' v ' 



d'^^q-qg ^ fcf(£l -£2) - 2(fct ■£i)(fct -£2) + 3qQ^ 



Moreover, since the mass M is formed from two subsystems with fc_|_ components ag+ 
and (1 — we have 



and hence 



M2 ^ {k^^k^f = {q+, kll{a{l~a)q+), (32) 



ki = Ar a{l-a). (33) 



Therefore for both the 'hard' and 'soft' regimes, in the limit of small diffractive masses, 
M <C Q, and keeping a fixed, we have for 

, 9^*^^ ^ kt ^ M, (34) 

^qg : ^^*gg ^ COUStaut, (35) 
^qqg , dHqqg ^ kl ~ (36) 

One also finds from the analysis of the wave function that the longitudinal part is a 
higher twist contribution, i.e. the longitudinal structure function is suppressed by an 
extra power in at fixed /?, where 

^ ^ ^ = OT- '''' 



One can show this result by substituting a in ( pO[ ) and by fc^/M^ according to 
(p3|) assuming that a is small. One then obtains for (pG) an extra factor kt/Q as 



compared to (29) 



X To be precise we should use the 'skewed' gluon distribution since the momentum flowing along the 
exchanged gluon lines in Fig. 3 are not exactly equal and opposite, see Section 11. 
f In the 'soft' limit the two gluons couple to the same quark line, that is they act like a single pseudo 
exchange particle with C = 1 and therefore give the same result as in (jl^ psj . 
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The photon wave function determines the general structure of the (3 spectrum. 
We have aheady noted that at large masses M the contribution to the inclusive 
difFractive process 7*p — > Xp arising from qqg production dominates over qq even 
though it is higher order in ag, since it contains gluon exchange as opposed to quark 
exchange (see Fig. 3). In summary, (p^)"(^6|) indicate that the /3-distribution has a 
distinct separation into three regions of small, medium and large (3 (large, medium 
and small M) where each of the three contributions, qqg, transverse and longitudinal 
qq respectively dominates. The fact that the longitudinal part is non-negligible also 
indicates the importance of higher twist terms even at fairly large . One of the 
important conclusions in Ref. |2J] was the observation that the wave function for qqg 
results in a rather soft gluon distribution (strongly decreasing as /3 1), see (|36|). 
It has been shown that a parameterization based on the wave function formalism 
provides a good description of the data without the need of a very hard, singular 
gluon distribution (strongly peaked as /3 ^ 1) as suggested by HI p6| . 

We use the it factorization formula ( p8|) to obtain the explicit form of the diffractive 
structure functions. First we note that the gluon distribution T is independent of the 
azimuthal angle of ^t, as t ^ 0. We can therefore easily perform the integration over 
this angle. We then take the square of the amplitude and integrate over t assuming a 
simple exponential form, exp(— where the diffractive slope Bjj is known from 



experiment. The final result for the diffractive structure function F2 
the sum of the three contributions 
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We have changed variables in the square brackets using ( |3^ ) , and in addition we have 
introduced 



e = 



J±_ 

l-(3 

l-z 



for qq, 
for qqg 



(41) 
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Note that the infrared cut-ofF is hidden in the unintegrated gluon T of (p^). It is 
specified by the size of the proton. Eqs. (|3^-|40|) are therefore safe in the infrared 
Umit, fc^ ^ 0. The variable z appears for qqg instead of /3 because we stih have to 
convolute with the quark box. The convolution is apparent in terms of the Altarelli- 
Parisi splitting function, describing the g qq transition (that is the second factor in 
the second line of (|40|)). The z describes the relative momentum fraction of the gluon 
with respect to the Pomeron. 

The unintegrated gluon distribution is the only quantity in ( p8[ - ^ ) still to 
be determined in order to perform a numerical analysis. It can be obtained from 
standard parameterizations for parton distributions or by using a particular model as 
in Ref. [p2[ . Fig. ^ shows how the three contributions, qqg, qq from transversely and 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 



P P 

Figure 4. The diffractive structure function and its three components plotted 

versus (3 for a fixed xp = 0.005. Dashed line: transverse qq, dashed-dotted line: 
longitudinal qq and dotted line: qqg. 



longitudinally polarized photons, occupy the three different regimes at low, medium 
and high /3. This decomposition is in large due to the property of the wave fuctions 
as argued previously. The precise form of the parameterization for J- has only little 
influence on the P spectrum. The relative normalization is given by colour factors. 
Since the colour factor for qqg is much bigger than for qq, it is another reason why it 
becomes important despite the fact that it is a higher order contribution (as was set 
to 0.2). One notices the drop of the curve for at large (3 when is increased. This 
effect is due to the higher twist nature of the longitudinal contribution as discussed 
earlier. The rise at small f3, on the other hand, is driven by the logarithm in which 
results from the phase space integration of the quark box. 

The X ip-distribution depends very much on the parameterization that one chooses 
for the gluon distribution In general, since we are considering a perturbative 
approach, the slope is much steeper as compared to a soft approach. As an example for 
a soft approach one can take the model of Ref. which is based on nonperturbative 
gluon exchange. It leads to a much shallower xjp-distribution [§8|, pTf. 
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5. Impact parameter representation 

It is informative to recast the formulae derived so far in impact parameter space. In 
this way we can see the physical picture described earlier in which the position of the 
partons in impact parameter space is frozen during the scattering. This property, and 
the fact that the colour-dipole approach to small x physics in Ref. 1 28 , 2^, ^ ^ lives 



in impact parameter space, has led to its increased popularity. Moreover, the concept 
of the dipole cross section can be generalized from two-gluon exchange to multi-gluon 
exchange. The disadvantage, however, is the need to transform back to momentum 
space when exclusive distributions such as the /3-spectrum are studied. 



We transform the wave functions ^'(fcf , at) of (p4[)-(26) to impact parameter space 
using the Fourier transformation 

^{r.a)= (42) 

where r denotes the separation of the dipole in impact parameter space. Carrying out 
the integration we find 

V'T-(r,a) = a^(l - a)^ Q K^{^/ a{\ - a)QV2) f • e(A^ = 1), (43) 

ZTT 

V'-+(r-,a)= - ^^a^(l-a)igifi(y^^(r^^^0Q^)^--£(A7 = l), 

ZTT 

7/.^_(r,a) =^^+(r,a) = ^ a(l - a) Q Kq{^J a{\ - a)QV2), (44) 



TT 

V'g^g(r, a)^ - 7^ (£i • £2 - 2r • £i r • £2) v^oQ^ K2{^J otQ'^r'^) , (45) 

ZTT 



where the Kn are modified Bessel functions of the second kind. We can use these wave 
functions to evaluate the diffractive structure functions. For example we find 

^ 64^55^ 4- Ml-/?)' J (27r)2 * Vl - 4a2fc2/Q' 



1 



(fr (fr' e''^'-^'^-'^'^ a{r,X]p)air',xip)r-r Ki{aktr)Ki{aktr'), (46) 



where = /3/(l — /?) and where we have introduced the dipole cross section 



47r 



air,x)^^J^ 1-e'^-^* asT{x,ei) 
47r2 /■ d£ 



3 y [1- Jo(M] asH^^f,) . (47) 

Inserting ( ^ ) into ( ^6| ) and computing the integral over r or r' one directly gets 
back to the expression in (^8|). a-(r, x) is the efi^ective inclusive cross section for the 
scattering of a qq system, with g, q transverse separation r, on the proton at energy 
y/s, where s — M'^/x. The diffractive or elastic scattering amplitude is given in terms 
of (T, via the optical theorem. Hence the aa structure of the observable (^). 

One important motivation for using the impact parameter representation is the 
factorization into the dipole wave function and dipole cross section, which is equivalent 
to fcf-factorization in momentum space. The impact parameter does not change in the 
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course of the scattering, i.e. the incoming and outgoing states have the same impact 
parameter. This feature is clearly visible in the /3- or M-integrated version of (^). 
When f3 is changed back to a one obtains 



da [a^ + (1 - af 



^T,,gl^'t^ )- i287t^Bd ^ « Jo 

X a{l ~a)Q^ J Kl{^a{l - a)Q'^r^) a^{r,x) (48) 

which involves the square of the dipole cross section a. That is r = r', and the 
impact parameter is unchanged by the interaction with the proton. Mass eigenstates, 
however, are not eigenstates of the impact parameter (see (^)) p2| . 



6. Factorization in j*p Xp and diffractive parton distributions 

A general proof of coUinear factorization in diffractive deep inelastic scattering has 
been given in Ref. Q (see also @). In Ref. it has also been stated that 
factorization is violated in diffractive hadron-hadron scattering. We will demonstrate 
this in section ^ In general we cannot factorize diffractive cross sections into 
a convolution of a 'hard' partonic subprocess with universal diffractive parton 
distributions. On the other hand the factorization proof allows us to introduce such 
distributions for diffractive deep inelastic scattering processes. 

We can demonstrate from the expressions for the diffractive structure functions, 
that we can extract diffractive parton distributions which are in agreement with 
collinear factorization. For qq production we only have to strip off terms from F-Tqq 
which are of higher twist nature (such as a^k^/Q"^ in (^6|)). This of course means 
that the longitudinal contribution is obsolete in this context. In direct analogy to the 
parton model of ordinary DIS, (^), for diffractive DIS processes we write 

= E^^/5,^(/3) (49) 

9 

where we have introduced the diffractive quark distribution. Notice that /3 plays the 
role of the Bjorken x variable for diffractive processes. Identification (^ ) is made for 
fixed values of Xjp and Q^. It is often said that q^{(3) is the quark distribution of 
the 'hard' Pomcron with the quark carrying a fraction /3 — x/xp of its momentum. 
Of course to introduce such a picture we would have to specify the probability of 
finding the Pomcron in the proton, and to assume that the Pomeron is a real particle 
(that is a pole in the t plane or complex angular momentum plane). The ambiguities 
and difhcultics of interpreting q^{/3) as the quark distribution of the Pomeron are 



discussed in Section 14. By comparing (49) with the leading twist part of (^6|) we may 



introduce the diffractive quark distribution (see also [ p5| ) 

_ ^ ^ f '^^^t 1.4 f ,2 f ,2 / tkfir~r') 

- USn^xpBn [I - J {2^f K ]dr jdr e 

X (T{r,xp) (j{r' ,xp) r ■ r Ki{aktr) Ki{aktr'), (50) 
where recall = /3/(l — /3). Similarly using the qqg structure function F^^g we may 
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introduce the difFractive gluon distributionf 

fjD(fi\ - ^ f z,4 [ ,2 f f t ikf(r-r') 

X a(r, xp) a{r', xp) (<5™" - 2f"f") ((5™" - 2f""f'") K2{aktr) K2{aktr') (51) 

with m,n — I, 2. There are differences between using the full expression ( p6[ ) for xpF^ 
and using of (pG) in (ES). One difference is in the form of the fcj integration; in 



particular in the absence of an upper limit on fcf in (50). This means that energy- 
momentum conservation is obviously violated. The integral is nonetheless well defined. 
Imposing the true kincmatical cut-off would 'only' lead to higher twist corrections 
which in the context of coUinear factorization are subleading. The importance of the 
longitudinal contribution in comparison with the data, however, already indicates 
that higher twist corrections are not completely negligible. Another issue is the 
factorization scale ff . A dependence of the diffractive parton distributions on the scale 

is introduced by evolution, where ( |50| ) and (jsj) serve as the initial distributions. 
The initial scale Qq, however, cannot exactly be determined, which introduces a 
considerable uncertainty for the prediction. 

In the conventional parton picture there is always a soft remnant which carries the 
opposite colour to that of the elastically scattered parton. For the dominant diffraction 
contribution the picture is the same. There are, however, events with no soft remnant. 
One example is the exclusive qq final state with large k^. Although subleading, these 
events constitutes an example which breaks factorization. An important signal for 
these events is the azimuthal angle distribution of the q and g jets 

This discussion is not intended to disprove factorization but to point out some of 
the limitations. Diffraction has the unique property that the whole event is contained 
in the detector, i.e. the 'soft' remnant is included. Factorization focuses on the hard 
subprocess and basically disregards the soft remnant. The approach we are pursuing 
based on Feynman diagrams can overcome part of these caveats. It offers a complete 
description of the final state where the remnant is represented by a quark or gluon. 

7. Triple Regge limit 

The qqg contribution in the triple Regge limit, i.e. ^ Q^, can be calculated 

assuming the strong ordering of the longitudinal momentum components. We can 
relax the ordering condition on the transverse momenta instead. A further requirement 
in this approach is the strong hierarchy between the large diffractive mass M and the 
transverse momenta kt of the quarks or the gluon, i.e. M ^ kt- 
The diffractive structure function in this limit reads [p9|: 



^^Fm = T.4 512^3^1 I' '^^ + - / ^'^^^'^2 ^^^^^^^^ 



(52) 



where i,j = 1,2 are vector components in the transverse plane. The amplitude is 
given by 

(f£ ^, „2. ( f £ + ki+k2 ki+k2 



r (p£ ( 



D{i + ki + k2) D{ki + k2) 



t The Operator Product Expansion allows diffractive parton distributions to bji_^introduccd 
consistently to any order in logQ^. According to ]36| the initial distributions of ( pcj ) and 
are valid beyond leading logQ^. 




Figure 5. Two typical diagrams contributing to diffractive qqg production. In the 
first diagram the momentum of the off-shell quark is ki + fc2, and in the second is 
fci + fc2 + £. The dotted line denotes a cut through the diagram (giving twice the 
imaginary part of the amplitude, which dominates over the real part). 

and where ki and k2 are the transverse momenta of the quark and the gluon 
respectively (see Fig. P). The first factor (• • ■)i in ( p3|) represents the different forms of 
the propagator of the off-shell quark in diagrams such as those in Fig. 5. The second 
factor {■ ■ ■)j is associated with the coupling of the s channel gluon (to be precise it is 
the BFKL vertex in thep-A = gauge, which also encompasses gluon Bremsstrahlung 
emission); it also includes the gluon propagator. Unlike (p6|), here a is the longitudinal 
momentum fraction of the quark and not the gluon. The longitudinal contribution to 
has been neglected here (see Ref. [|9j). One can also study the dependences on 
the azimuthal angle of fci and k2. 

For practical purposes the representation in momentum space is inevitable. The 
previous result becomes more compact, though, if one transforms ( |5^ ) into impact 
parameter space: 

J (fki(fk2 M.,jM*j = 324 a(l - a)Q^ j d?rid?r2 (^^J a{l - a)Q^rl 

^ \^2 ['^(^2) + Hr2 - ri) - a{r,)f . (55) 

ri {r2 - ViY 

The vector ri represents the separation of the quark and antiquark, r2 the separation 
of the quark and gluon and r2 — Ti the separation of the antiquark and gluon. The 
combination in which the dipole cross section appears in (p5|), 

'7(r2) + <7(r2 - ri) - <5-(ri), (56) 

refiects the scattering of the three effective colour dipoles, qg, qg and qq, on the 
target proton. Note that the scattering of the gg-dipole in this configuration in not 
colour suppressed by powers of N,.. The reason for this is a contribution in which 
the initial gg-pair interacts before the gluon is emitted (second diagram of Fig. 
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Similar results have been found in Ref. p30| . In the short distance limit, ri ^ r2, the 
qq-paiv recombines effectively into a gluon, forming one pole of the gluon-dipole. The 
scattering is given by a single cross section, (j(t'2), for the gluon-dipole. 

The result for the triple Regge limit is consistent with previous results when, in 
addition to the ordering of the longitudinal components, the ordering of the transverse 
momenta is assumed. In this limit (/3 —i- 0) one can show that ( pO[ ) and ( p^ ) coincide. 

We have originally argued that the inclusion of a qqg component of the photon 
wave function at large M is enough for most purposes. In the limit of really high 
energies, however, with each order in as one gains a logarithm in M'^/Q'^ which can 
overcome the smallness of the coupling. The resummation of the ctg log(M^/Q^) 
terms corresponds to the emission of multiple s channel gluons. In Ref. [|40[ it has 
been shown that the t channel structure then becomes complicated. The two two- 
gluon exchange ladders in the diagram for the diffractive structure function interact 
leading to a 4-gluon "bound state" in the t channel. There is no simple coupling of 
three gluon ladders with a local triple 'BFKL' ladder coupling which would enable 
the two exchange ladders to coalesce into a single ladder. Nonetheless, one can find 
special kinematical configurations such that a triple ladder configuration exists ||43|| . 
This happens at t = when the upper ladder is forced into the DGLAP regime 
where leading logs in become dominant. The lower ladders on the other hand 
are of 'BFKL' type. This result is due to a certain scaling behaviour which leads 
to the 'conservation' of anomalous dimensions of the three ladders. For i ^ 0, 
however, this scale behaviour is absent and the coupling of three BFKL ladders is 
in principle possible except that the 4-gluon bound state supersedes the simple triple 
ladder scenario. 

8. Multiple gluon exchange and the semiclassical approach 




Figure 6. Reduction of multiple- gluon exchange into eikonal factors shown by the 
dashed lines. The colour structure is discussed separately. 



So far as we have discussed diffractive processes using the perturbative assumption 
of two-gluon exchange. Is perturbative QCD really applicable, particularly as it 
embraces potentially 'soft' contributions where as becomes large? In the soft regime 
we would expect multiple gluon exchange. In QED the generalization from two- to 
multi-gluon exchange is fairly straightforward. One simply has to introduce eikonal 
factors 1 17 which account for the effect of resumming multiple photon exchanges, 
while the photon wave function does not change. One can visualize this procedure by 
merging any number of photons coupled to one arm of the dipole into one effective 

t Loosely speaking, BFKL cvo ution is the resummation of log(l/a;) terms applicable at small 
X, whereas DGLAP evolution corresponds to the resummation of aglogQ^ terms. 
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vertex (Fig. In QCD, due to non-commutation, the colour structure turns out 
to be non-trivial. However this can be absorbed into the parameterization of the 
unintegrated gluon distribution. What matters is the relative colour factor of the 
quark dipole compared to the gluon dipole (see Fig. 3). Taking the limit of a large 
number Nc of colours, this is the same as for two-gluon exchange. The argument is 
that the leading- A^c colour tensors, for a gluon loop with an arbitrary number of gluons 
attached to it, can be decomposed into a set of tensors with coefficients depending on 
Nc [Q. The leading tensor is identical to the quark loop tensor times Nc- The quark 
and 'gluon' loops simply come from squaring the qq and qqg production amplitudes 
of Fig. 3. We have illustrated this for four gluon exchange in the first row of Fig. 0. 
In this scenario one indeed finds that multiple gluon exchange couples to a gluon or 
a quark dipole very much like a simple two-gluon exchange, i.e. all the parametrized 
gluonic structure of the 'Pomeron' is embodied in the distribution T. 




Figure 7. Relating the colour structure of the cross section for {qq)g production to 
that for qq production in the large Nc limit for the example of four- gluon exchange. 
In the Glauber approach the factor is 2^ Nc where P is the number of ladders 
exchanged, whereas in the general case the factor is always Nc in the leading Nc 
limit. 

We may confront the general approach with a frequently used simplified model 
in which it is assumed that the gluon ladders are non-interacting. This model is a 
specific subset of the general t channel gluon structure and it is often assumed, though 
as yet unproven, that this represents the dominant configuration. It is associated 
with shadowing and unitarity corrections and so we call it the Glauber model. In 
this approach where the exchange of multiple colourless gluon pairs is resummed 
the situation looks different to the, more general, interacting ladder scenario. The 
projection on pairs of colour singlet states mixes the leading Nc tensor with non- leading 
tensors. With regard to the overall counting of powers of asNc the Glauber approach 
is subleading. However, the only effect of summing the leading asNc terms (associated 
with the leading- A^c colour tensors mentioned earlier) might be the reggeization of the 
initial gluon pair and would, therefore, not contribute to unitarization or shadowing. 
In this case the Glauber type multiple scattering would gain importance. As the 
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second row of Fig. |^ shows, once the four gluons are projected pairwise on a colour 
singlet state the relative colour factor of a gluon loop compared to a quark loop is 
ANc instead of Nc- The factor 4 is the large Nc limit of the ratio C\/Cp, i.e. for each 
pair of gluons we find a relative factor Ca/Cf for the scattering of a gluon dipole as 
compared to a quark dipole. In this scenario the parameterization for the unintegrated 
gluon distribution J- would be different for qq and qqg production. Both contributions 
would no longer be related as in (|3|) and (Eo|). 

The semiclassical approach in Refs. |35|, p5[ is conceived as a non-perturbative 
approach (for a recent review see Ref. |Q) and utilizes Wilson loops. It is very much 
reminiscent of the eikonal approach in QED. According to our preceding arguments 
the photon wave functions turn out to be the same as in the two-gluon exchange. The 
differences emerge in modeling the dipole cross section. The semiclassical model shows 
similarities to the Glauber approach and has therefore not quite the structure of ( ^8|) 
and (^0|). The difference occurs in the treatment of the colour as we have discussed 
above. In contrast to the Glauber approach the energy dependence is introduced as 
part of the overall normalization, which is taken to have the form ^1{L — \n x)^ , where 
and L are parameters to be determined by the data. It is argued that this form 
of energy dependence is due to the number of field modes which increases when the 
energy is increased. It is consistent with unitarity constraints. 



9. Colour dipole approach 

The colour dipole formalism has been developed in |2^, ^ as an alternative 
to the Feynman diagram approach to small x physics. It is formulated in impact 
parameter space and has been shown to reproduce Feynman diagram results for 
inclusive processes in the Regge limit, as embodied in the BFKL equation (which 
will be introduced in Section 13). With regard to gluon radiation in diffraction it 
can be applied in the triple Regge limit, i.e. for large masses 3> only. At 
lowest order and first non-leading order in the strong coupling, i.e. qq and qqg final 
states, the colour dipole approach is consistent with the results derived in Section 7. 
In particular it can be proven to coincide with the triple Regge result of (55) 



One expects that multiple gluon radiation in the colour dipole approach leads to the 
same results as the Feynman diagram approach po| . A complete phenomenological 
treatment of diffraction within the framework of the colour dipole approach has been 
presented in Ref. |Q . The colour dipole approach of Ref. jis) uses the non- forward 
BFKL-Pomeron and, in this respect, goes beyond other approaches. It directly delivers 
the t integrated cross section. 



10. Diffractive production of open charm in DIS 

The measurement of charm in diffractive scattering provides an additional test for 
any approach based on the photon wave function formalism. As before in the case 
of light flavours we consider the exclusive cc-pair which arises from the dissociation 
of longitudinally and transversely polarized photons, as well as the production of 
the ccg-state. Since the mass of the charm quark sets a limit on the size of the cc- 
dipole, it becomes colour transparent and one expects a strong suppression for this 
configuration. The effective gluon dipole, associated with ccg production, on the other 
hand, is not restricted in size. 

For a quantitative analysis we have to extend our expressions for qq production 
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to include the charm mass m^. 
7*p ^ ccp are M, g 



We find that the difFractive structure functions for 
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The definition of /3 is the same as before, i.e. /3 — Q^/{M^ + Q^), but the diffractive 
mass M has now a lower limit given by 2mc. Hence the P spectrum has an upper 
limit well below 1 at the lower values of Q^. Also the two-body kinematical relation 
( ^ ) acquires a mass term and changes into 

2 + Trie /I \ fai \ 

^=^^-"(i-«)y (61) 

Note the second term in the curly brackets of (^7|) which results from a spin flip on 
the quark line. The spin flip is only present in the massive case. 

We may go a step further and investigate the limit of small £f in the expressions 
( p9| ) and (|60|). The relevant scale for the expansion is k"^, i.e. we expand in powers of 
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(62) 



(63) 



Here g is the conventional gluon distribution of the proton, obtained by integrating 
the unintegrated distribution The interesting point about the scale k"^ is its lower 
cut-off m^/(l — (i) as kt approaches zero. This means that the diffractive production 
of exclusive cc-pairs always stays in the hard or perturbative regime regardless of 
the value for kt. At large fc^, i.e. k1 ^ mj,, the distribution in k^ falls roughly as 
1/kt for the transverse part. We can therefore approximate the /cf -integration for the 
transverse part by replacing the argument k'^ in the structure function by mj^/ (1 — /3) 
and perform the integral for the remaining terms (see also Ref. [^): 
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(64) 



The cross section shows a rapid growth proportional to [xpg^ as long as the 
saturation regime is not reached. The prefactor l/?7i^, on the other hand, leads to 
a suppression as compared to light flavours (colour transparency). The longitudinal 
part has a /c^-spectrum which is proportional to l/k^ so that a logarithm builds up 
after integrating over k^■. 
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(65) 



It has already been noted that the longitudinal part is a genuine hard contribution 
which can be calculated using ( ]30| ) . The price to pay is the factor 1 /Q^ which makes it 
a higher twist contribution. The dependence of the unintegrated gluon distribution J- 
on fcj cannot be completely neglected. In this sense the previous equation represents 
only a crude estimate. 

For the ccg component we will proceed in a slightly different way, making use of 
the diffractive factorization property. In (|5l| ) we have specified the diffractive gluon 
distribution which needs to be folded with the corresponding charm-coefficient function 
191,0: 
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where a = 1 + 4m^/Q^ and 
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V, the centre-of-mass velocity of the charm quark or antiquark, is given by 
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In Ref. ||50|] the diffractive gluon distribution was calculated assuming a cut-off 
fco on kf (see also Ref. |5^). This result can be rederived from the momentum 
representation^ of g^ by expanding in powers of it/k'^: 
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bearing in mind that fc^ — k^/{l ~ (3). Another way of deriving the second equahty 
would be the use of the second derivative (^l|) of the '^qqg wave function. Unfortunately 
only the hard regime is being correctly taken into account. More appropriate is the 
use of the full expression in the first line of (^9|) . 

Another comment is appropriate, which concerns the possible K factors suggested 
in Ref. [Q. K factors are used to represent the higher order corrections to the 
process. A certain virtual correction has been singled out which in an extended 
analogy to inclusive Drell-Yan pair production might lead to a significant enhancement 
of the prediction for diffractive scattering in general, not only for diffractive charm 
production. In order to understand the importance of higher order corrections a 
complete NLO calculation is required. 

The particular signature for diffractive charm production is a fairly large 
contribution around 25% at small /3 where the ccg component dominates. The fraction 
of charm in this regime is the same as expected in inclusive charm production. At 
larger (3 and large Q^, i.e. /3 > 0.2 and > 10 GeV^, however, charm is suppressed 
due to colour transparency. It contributes only a few percent, with a possible K 
factor enhancement. This feature helps to discriminate against other approaches (for 
example in Ref. pof) which suggest a hard and dominant gluon distribution at large 
/3 and consequently a much larger charm fraction. 

11. Diffractive vector meson production 

The diffractive electroproduction of vector mesons, "f*p Vp offers the opportunity 
to study many aspects of perturbative QCD, as well as revealing information on the 
gluonic structure of the proton. The hard scale can be either the virtuality of 
the photon the mass of the quarks for heavy vector mesons, the momentum 
transfer t or some combination. Data are available for p, w, (j), J/ip, ip' , Tf" vector meson 
production. Moreover the observed decays such as p ^ tt+tt" and J/V' M^p~ 
allow a helicity decomposition of the 7*p — s- Vp amplitudes. The basic mechanism 
is the same as that for open qq production, except that we must now include a wave 
function tp^ {qq) for the eventual formation of the vector meson, see Fig. 8. Since, 
at high energies, the 7* — > qq fluctuation time and the qq V formation time are 
both much longer than the qq interaction time with the proton, the amplitude has the 
factorized structure 

A{-f*p ^ Vp) = ® Aqq+p ® 4,^^. (70) 

We must also note that 7*p Vp is described by "skewed" parton distributions 
of the proton, that is x,lt 7^ x' ,1^ in Fig. 8. Reviews^ of "skewed" (also called off- 
forward, non- forward or off-diagonal) distributions can be found in [ p7| , p8| , |6^. For 
t = the skewedness comes from x ^ x' , and increases as the mass of the vector meson 
Mv or the photon virtuality increase. However, for small x, relevant to 7*p Vp at 
high energies, the skewed distributions are completely determined by the conventional 
diagonal distributions |6l| . Since the cross sections depends on the square of the gluon 
distribution, these processes are expected to be a particularly sensitive constraint on 
the gluon. The subtlety is the determination of the vector meson wave function. 
Various approaches give quite different results. The corresponding uncertainty mainly 



t There is much recent activity, and a large literature, concerning skewed parton distributions. For 
the high energy diffractive processes that we discuss, we are interested in distributions in the small x 
domain. A general review focusing on their use in virtual Compton scattering can be found in |p6|. 
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affects the absolute normalization, and less the energy dependence of the cross section. 
But even the energy dependence is not perfectly well established, for the scale which 
enters the structure function is not linked to the hard scale in a simple way. 




V = p or J/ip 



Figure 8. One of the four basic diagrams for diffractive vector meson 
electroproduction, 7*p Vp. (There are three other couplings of the gluons to 
the qq pair.) In general x ^ x' , It ^ t'^ and the cross section depends on the square 
of the "skewed" gluon distribution. 



On the other hand for large t the skewedness in the gluon ladder comes from 
£t 7^ £[ and the evolution follows BFKL. Large t diffraction is an ideal place to study 
the high energy limit of pQCD Unfortunately large t vector meson production is 
again hampered with the problem of how well the meson wave function is understood 
psf . It has therefore been proposed to look for photons instead of vector mesons in 
the final state ||2^, i.e. quasi elastic 7p-scattering. At t = 0, but large Q^, this 
would be Deeply Virtual Compton Scattering (DVCS) 

Vector meson production has developed into a very active field which deserves a 
separate review. Here we will concentrate on the perturbative QCD description of 
these processes, but note that there are also approaches which are based on the vector 
meson dominance (VMD) and related models |6^, First we will study J/ip and p 
production at t = and then, in the following section, t ^ production. 



12. Diffractive J/t/j production at t = 

We start our discussion with the electroproduction of J ftp |68|, p9|l. The modeling of 
the J/ip wave function is in general a non-trivial procedure |7C|, [72| but, for the 
purposes of illustration, we use the nonrelativistic, static wave function assumption 
made in Refs. |69| . 

The production mechanism for vector mesons is closely related to open diffraction. 
The photon dissociates into a gg-pair which then scatters off the target. With a 
certain probability the final state recombincs into a vector meson. This probability is 
determined by the projection on the meson wave function. The photon dissociation 
is as before described by the photon wave functions ( p^ - |25| ). We can make a 
simple qualitative estimate about the longitudinal versus the transverse cross section 
assuming that is much bigger than {M^ is the mass of the J/i/') and the 
corresponding [3 is close to 1. Whereas the transverse component^is decreasing when /3 

f Note that the expression in the curly brackets in (^) vanishes as /3 — > 1, and overcomes the 
prefactor 1/(1-/3). 
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is approaching 1 the longitudinal component remains constant. This already indicates 
the dominance of the longitudinal versus the transverse cross section. 

The wave function of a vector meson is constructed in analogy to the photon wave 
function. In the simplified nonrelativistic, static case one assumes that the quark and 
antiquark are weakly bound in comparison to their mass and that, therefore, their 
relative movement is small. This means that the relative transverse momentum is 
approximately zero, i.e. /ct « and that the quark and antiquark equally share the 
longitudinal momentum fraction, i.e. a ~ 1/2. It also means that the mass of the J/ip 
is approximately the sum of the mass of charm and anticharm quark, i.e. = 2mc. 
In this simplified approach one can write 

*^(fct,a) - 6ia~l/2)d^ikt) . (71) 

We take the expressions for open charm production, (^) and (58), and apply the 
(5-functions of (^l]). To this end one has to integrate over (3 using the relation ( |6l|) 
which connects a with /?. The first thing to notice when kt is set to zero is that the 
integral It of (|59| ) vanishes. The final state configuration has to match the spin of the 
incoming photon. For a transverse polarized photon this means that the final state 
has to have the total angular momentum = ±1. This can be the 'orbital' angular 
momentum of the emerging charm- anticharm pair, — ±1, when the quark spins are 
opposite, i.e. Sz — 0. The other possibility is spin-flip along the quark line, 5*2 = ±1, 
while the orbital angular momentum = 0. For the orbital angular momentum to 
be non-zero the transverse momentum kt has to be non-zero. Since we require kt = 
we enforce Lz — and It = 0, leading to the vanishing of the first term in the curly 
bracket of (|5^). The second term, which is associated with the spin- flip, remains and 
leads, after a httle algebra, to: 



{Q' + M'^r \j * et+Q' 

IQQ^Ml _ 2 

- (Q2^M|)5 [^sxpg{xp,Q )\ (72) 

where xp = (Q^ + M^,)/{Q^ + W^), see (|) and (§, and 

= a(l - a)Q^ +ml ~ {Q^ + Ml) /A. (73) 

The leading contribution to the (."l integration comes from the region t1 < Q^. 

For longitudinal polarized photons both the orbital angular momentum and the 
spin are zero, i.e. Lz — and Sz = (no spin- flip). There is no suppression when kt 
is set to zero 

,2 asT{xp,ft] 



{Q 



l^^M^)^ ['^sxpg{xp,Q'')]' . (74) 



We see that, in this simple approach, the transverse and the longitudinal production of 
J /ip mesons are identical up to a factor Q'^ /M^ by which the longitudinal is enhanced 
over the transverse component. Also the dependence on the square of the gluon 
distribution is evident. 
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With regard to the absolute normahzation one has to go further into the details 
of the construction of the vector meson wave function. For reasons of consistency the 
same wave function determines the leptonic decay width for the J/ij} |^8[ |7^. This 
allows the normalization to be determined by the decay width T"^^ _ . The complete 



expression for the cross section (longitudinal and transverse) then reads 69, [72| 



The relativistic corrections to the J/V' wave function (|7l|) have been considered in 



Refs. [71 
Refs. [72 



72, 73 . The change that these make to the cross section is under dispute. 
73 [ show that the net effect of the corrections is small. On the other hand 
the authors of Rcf. [|7^ argue that they cause a large suppression. These corrections, 
together with the estimates of the higher order contributions, lead to an uncertainty 
in the absolute normalization||, but have much less impact on the predictions for the 
dependence of the cross section on the 7*p energy W . 

In the photoproduction limit we simply set = Q. Only the transverse component 
survives in this case. The large charm mass provides a hard enough scale for 
perturbation theory to be applicable. The experimental advantage of photoproduction 
is the high statistics compared to electroproduction. The analysis of Ref. [Q 
demonstrates the value of J /ip photoproduction in determining the gluon distribution 
of the proton. Although corrections which go beyond leading log in have an 
influence on the cross section, these mainly affect the absolute normalization. The 
shape of the W distribution seems invariant under these corrections. The prediction 
is 



\x 



using the rise of the gluon distribution with decreasing x determined by analysis of all 
available deep inelastic scattering data. This should be compared with the 'soft' or 
Pomcron Regge exchange prediction 



a 



where the Pomeron trajectory ap is given by ( [T^ ) and (• • •) indicate that the 7*p — > Vp 
cross section is averaged over t. We see, from Fig. ^, that the measured energy 
dependence of J/V' production is in good agreement with (H), which provides support 
for the perturbative QCD approach. 



13. DifFractive p electroproduction at < = 

The diffractive production of p mesons, 7*p — > pp, is special because this is the process 
for which the data are most extensive at present. However, unlike J /ip production, we 
do not have the mass of the charm quark to provide a hard scale. The light u and d 
flavours do not provide the hard scale and so only electroproduction can be considered 
for a perturbative QCD treatment, where large provides the scale. In fact the 
7*p pp data in Fig. ^ show the transition from the 'soft' Regge prediction of ( [77| ) 

f The use of a skewed gluon distribution causes a 30% enhancement of the cross section for J/i/i 
photoproduction [pl|. 
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Figure 9. T/ie energy dependence of the cross sections for diffractive p 
electroproduction (at three values of Q^) and J/ip photoproduction. Also shown 
is the total 7p cross section compared to the iy0-l6 energy dependence of (u4j- The 
W^'^^ and ly"'* lines indicate the 'soft' Regge Pomeron and 'hard' QCD Pomeron 
predictions of and ( [t^ respectively. (A. Levy, private communication.) 



at Q2 = to the 'hard' QCD prediction of (|76|) at w 20 GeV^. The other evidence 
that perturbation theory is apphcable at large comes from the measurements of the 
slope Bp which defines the diffractive peak, da/dt ~ exp{—Bp\t\). The value decreases 
rapidly from Bp - 10 GeV^ at = to Bp - 5-6 GeV^ for - 10 GeV^. The 
latter value is expected from the size of the proton which suggests that the size of the 
"f* p vertex is close to zero. 



From (75) we see that the leading order prediction for electroproduction in 



longitudinally polarized states is 

[xgix,Qy4)r 1 
^ 

for Q2 ^ M^, where x = Q'^ /W^ and 7 is the effective anomalous dimension of the 
gluon density 

xg{x,Q') ^ (Q2)7. (79) 



We have taken a representative value 7 = 0.3 determined from the global parton 
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analyses, corresponding to the x range (10 ^; a; < 10 ^) of the HERA data. The 
QCD prediction (vm for is consistent with the behaviour seen in the data. This 



is not the case for ut- The non-relativistic prediction for ctt, obtained from (72) and 

= ^ OL, (80) 

appears to be too smah and to faU too rapidly with increasing , as is evident by 
comparing the dashed line in Fig. ^ with the data. The problem lies with the wave 
function which describes the recombination of a qq pair into the p meson [f70| . Unlike 
J production, it is crucial to include relativistic, non-leading twist corrections even 
to obtain a rough estimate of (Tl/ct 0|- However these turn out to be insufficient 
to resolve the problem. Moreover there are indications that a non-perturbative 
approach will make ut fall off even faster with and make matters worse. 




Figure 10. The dependence of the QCD predictions J7q/ for the ratio o-Zj/ctt 
of the electroproduction of p mesons {-f*p — > pp) in longitudinal and transverse 
polarisation siatei compared with HERA data. The dashed line indicates the 
expectations of (iSClj- 

In order to circumvent the difficulties of modelling the p meson wave function an 
alternative approach has been proposed Q . It is based on the production of uu and dd 
pairs in a broad mass interval containing the p resonance. In this mass interval, phase 
space forces the qq pairs to hadronize dominantly into 27r states. Moreover provided 
the gg-proton interaction does not distort the spin, we expect that the transition 
J* —> qq ^ 2tt will dominantly produce 2tt systems of = 1~. Indeed this is found 
when the = 1~ state is projected out. This parton-hadron duality approach allows 
detailed predictions to be made for p electroproduction. 
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The predictions may be obtained by using the resuhs of (|57|)-(|60|) to calculate the 
production of uu and dd pairs of mass M via two gluon exchange. We put nic = 
and change the integration variable in ( |57| ) and ( |58| ) from fcj to cos 6, where 9 is the 
polar angle of the outgoing q in the qq rest frame with respect to the direction of the 
incoming proton. Thus the transverse momentum of the outgoing quark is 

kt = iM sine*. (81) 

Then ( p7| ) and (|5^) can be shown to give the qq production cross sections |7j] 

dar _ 47r2e^Q: l 
dM2 ~ 3Bp (Q2 + M2)2 4 j_i 



y'^dcos0(|dii(0)|' + |d}_i(f?)|') (83) 



where d'l^ {d) are the usual rotation matrices. The di^ integrals It.l over the transverse 
momenta ±£t of the exchanged gluons are given by (^) and (|6^) . They are essentially 
the cross sections for the interaction of the {qq)L.T pair with the proton. If the qq +p 
cross section were constant (that is 1^ = It = constant) then we can perform the 



integrations in (82) and ( |83| ) and find aL/^T ~ /■iM^. An idea of the effects of the 
distortion of the qq state by the two-gluon exchange can be obtained by evaluating It 
and II of ( |59| ) and ( |60| ) in the leading log approximation, in which it is assumed that 
the main contribution comes from the interval < k^. Then 

II ^ It ^ {kjV/kl ^ (sin^ 0)^-1, (84) 

where 7 is the anomalous dimension of the gluon, see (|79|). If we substitute this 
behaviour in ( p2| ) and (|8^), project out the spin 1 component and assume that 7 is 
constant over the region of integration then we obtain 



(85) 



CTT M2 \7+l^ 

Now higher means larger x and both changes imply smaller 7. As a consequence, 
the decrease of 7 with increasing has the effect of masking the Q'^ growth of fj^ /err. 
The — \~ projection integrals over cos 6* turn out to be linearly dependent on the 
Ij (i.e. not on |/j|2 as in (|8^ ) and (|8^)) and are therefore less infrared sensitive than 
( ^ and (|8^). or, as well as ctl, is convergent as 6* ^ provided only that 7 > 0. 

The complete calculation does not, of course, make the above simplifications. The 
latest predictions ||75|| , which are shown by the continuous curve in Fig. uses 
a skewed unintegrated gluon distribution, which is determined from a conventional 
gluon found in a global parton analysis. Soft gluon emissions and the contribution of 
the real part of the amplitude^ are also included, ctl is found to be insensitive to the 
infrared, whereas ut is more sensitive but the uncertainty is less than the variation 
due to using gluons of different global analyses. The absolute normalization, on the 
other hand, is again subject to uncertainties. One is the higher order corrections which 
give rise to a substantial K factor, the second is the size of the AAf interval taken 
about the p resonance. Fortunately these ambiguities do not affect the ratio aj^juT- 
Moreover the shapes of the and W distributions can be predicted reliably. The 

t The optical theorem determines the imaginary part and a dispersion relation then provides an 
estimate of the real part. 
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enhancement effect of using tlie skewed distributioE||| increases witti Q^. Tlie QCD 
prediction for tlie slrape agrees weli with, the HERA data ||7^ . 

So far the discussions have assumed s channel heUcity conservation (SCHC) which 
means that the produced p meson retains the helicity of the incoming virtual photon. 
However a similar QCD model to that discussed above, predicts a small violation of 
SCHC given by the helicity-flip amplitude describing the 7t — > PL transition, which 



satisfies 78 



helicity-flip amplitude \/\t\/2 

(86) 



helicity-conserving amplitude Q 7 

The correlations of the observed angular distributions of 7* — > p production and 
the p — > 7r+7r~ decay allow 15 density matrix elements to be measured so that the 
consequences of ( p6| ) can be tested and, indeed, the fully helicity density matrix can 
be explored |7^. The QCD predictions Q are in excellent agreement with HERA 
data. 



14. Diffractive J/V' production at large \t\ and quasi-elastic scattering 

Diffractive J/tp production at large \t\ and even more so quasi-elastic 77? scattering 
are special cases of two particle elastic scattering as depicted in Fig. |ll| where the final 
states A' and B' might differ from the initial state particles A and B. The momentum 
transfer t is supposed to be large enough so that perturbation theory is applicable but 
still much smaller than the total energy W . As for inclusive and exclusive diffractive 
production at t = the driving mechanism is two gluon exchange at high energies 
where ^ \t\. Unlike the previous subsection where we had a two-scale problem 




Figure 11. Schematic representation of ^aa' '■^ ^f^^ impact factor for the 

A ^ A' transition via two gluon exchange, while G is the Green function for the 
two interacting (Reggeized) gluons. 



with a hard scale at the top (Q^) and a soft scale at the bottom we now have a 
single scale problem (t). In this domain contributions of the form as log(VF^/|i|) are 
dominant rather than as ^og{Q^ /Ql), and they must be resummed to all orders. This 
is accomplished by the non- forward BFKL equation |8^ . Let us sketch the application 

I An even larger enhancement, by a factor of about 2, comes from using skewed distributions to 
describe T photoproduction, 7p — > Tp 1761 [77|. 
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to elastic or difFractive scattering. The general structure of the amplitude for a process 
A + B ^ A' + B' is 

A{A + B A' + B') = <^AA' <»G(E)<i>BB', (87) 

as sketched in Fig. |ll|. The impact factors ^aa' and ^bb' describe the transition 
A A' and B B', whereas G is a universal Green function for two interacting 
gluons. The Green function, which we have previously called gluon ladder, satisfies the 
BFKL equation which has been derived and solved in the leading log approximation, 
where the leading (as ln(Vl^^/|t|))"' terms are resummed. In the next-to-leading log 
approximation as{as \n{W'^ / \t\)" terms would also have to be resummed. So far NLO 
corrections have only been performed for t — 

In the case of diffractive vector meson production A represents a photon (we will 
assume it to be real) and A' a vector meson. The problem of modelling the vector 
meson wave function is still one of the main obstacles. Since we are interested in the 
dynamcs of the 'hard' exchange rather than the final state, it has been suggested to 



replace the vector meson by a photon 1 20 , 8^ . Then A' is a photon again. The upper 
part of the diagram in Fig. ^ which describes the coupling of the two gluons to the 
photon has the same form as in ( ^0| ) except that we have to incorporate non-zero 
momentum transfer. The amplitude for elastic gg-proton scattering, A^g^p in (|70|), 
contains the Green function G and the lower impact factor In most cases the 

impact factor is defined in momentum space but for reasons which become clearer 
later on we work in impact parameter space. 

The amplitudes describing^ 7(7 — > ^jjq and 7(7 jq at large \t\ are of the form of 

(§) M 

A(7 + .^(^,7)+.) - j'^^J^^T^. {wh^T' ^''^ 

where My = 0, and Vs denotes the total energy at the parton level, since the 
exchange gluons couple to partons rather than to the proton as a whole, uj is the 
'famous' BFKL exponent |0 

Lo{v) = ^ [2^(1) - *(i + iJ/) - (89) 

where is the logarithmic derivative of the Euler gamma function, '^{z) = r'(z)/r(z). 
It is the key result from the resummation of the (as log(s/|i|))" terms by the BFKL 
equation. It controls the power behaviour with regard to s, that is the energy 
dependence of the process. This is why we speak of the QCD or 'hard' Pomeron. 



The general form for the 'impact' factor ^aa' for the upper blob of Fig. 11 is 



~ y"dVid2r.2Vnr-)S,(ri,r2) V'^^^*(r)e'"9-^i e'(i-")9-^^ (90) 

with r — ri — r2 and — \t\ and where il:{r) is given by (^2|). Since now \t\ ^ 
we have to introduce two impact parameter variables ri, r2, one for each of the quark 
and the antiquark produced in the 7 — *■ gg transition with longitudinal momentum 

f The reason why it is sufficient to consider the scattering off a quark, rather than a proton, is given 
above Q. 

X For pEoton production there are actually two different impact parameters <I>^-y. One comes from 
^'^{X'y = ±1)1/'''* (A7 = ±1) and one from the helicity flip product ip'^(Xy = ±1)1/'''* (A7 = =Fl)- 
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fractions a and I — a respectively. An integration over a is implied. $ is not quite 
the impact factor ^aa' of (p7|). First, part of the Green fmiction is included by the 
eigenfunction 

E^ir^,r,) = [£^y^^^ (91) 

of the non- forward BFKL equation and, second, it involves the convolution ((g)) over 
the impact parameter variables ri,r2. The effect of the net transverse momentum q 
flowing in the t channel is embedded in the phase factors in (90). The wave functions 
■0'''(r) and tjj'^{r) of the photon and J/tp are given, respectively, by ( |4^ ) and the 
'massive' version of (Q). Note that after taking the Fourier transform of "^'^{kt) of 
( [n] ) with respect to fct, we are left with S{a — 1/2). 

Large \t\ has the advantage that factorization works on the proton side. The 
Pomeron couples to individual partons rather than the whole proton. This statement 
is not as obvious as it seems for an exchange in the form of BFKL-ladder diagrams. 
It has been proven, though, in Refs. 0. The impact factor with respect to a 

quark-line reads: 

1.^,^2 Jd'r, (^^j e^1-^K (92) 

The factor 2 in front reflects the fact that either of the two gluons in the ladder 
diagram couple to the quark. For a gluon line the diff'erence is only the colour factor. 

Using the amplitudes ( p8| ) we may calculate da/dt for both jq — > i/jg and 7(7 jq 
at large \t\. The major work goes in evaluating the integrals in the impact factors ( pO| ) 
and (|9^). The numerical analysis of both cross sections shows a strong enhancement 
due to the BFKL resummation of leading (as log(s/|i|))" terms. The analysis also 
shows that both cross sections are of the same order of magnitude and very rapidly 
rising with increasing energy. The expected statistics at HERA should allow an 
experimental study of these processes. An important issue are NLO corrections to 



the BFKL-kernel |83|. For t = they were found to be very large and negative. 
For t 0, however, these corrections are not yet known. Nevertheless the general 
expectation is that they will reduce the theoretical prediction for the cross section. 
The quasi elastic photon-quark scattering can also be studied with virtual photons in 
the initial state. This would reintroduce the scale Q]. Once the impact factors are 
calculated they can be recombined in various ways. One example would be elastic 77- 
scattering at high energies. The limit t ^ can be studied in the case of 7*7* 77. 
Deeply Virtual Compton Scattering (DVCS) i.e. the process 7*p jp, is the 
analog process in deep inelastic scattering. The limit i — > as such is perturbatively 
defined, the t dependence is sensitive to non-perturbative effects. 



15. Non-factorization in diffraction at hadron colliders 



Historically diffractive scattering was measured in hadron collisions first. In the case 
of single inclusive diffraction one of the hadrons is tagged in the detector whereas the 
other dissociates. This process itself is, of course, purely soft, but with increasing 
energies the opportunity has emerged of measuring jets in the final state or similar 
hard processes (such as heavy flavour or VF-boson diffractive production). 

The first idea of how to interpret hard diffractive processes was to view the 
Pomeron as a quasi-real hadron with a parton structure that allows the description 
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of hard scattering in analogy to usual hadron-hadron scattering. One only needs 
to define parton distributions for the Pomeron, to be determined by experiment, 
and a Pomeron flux-factor. The latter represents the probability for a hadron to 
radiate off a Pomeron. This probability, of course, should be independent of the 
hard subprocess and consistent with soft inclusive scattering. One can assume the 
following factorization of the soft diffractive cross section into a flux factor fp and 
the hadron-Pomeron cross section cr^ : 

^fjp{xp.t)a^{M\ (93) 



dtdxp 

The variable xp has the same meaning as in previous sections (i.e. xp ~ h'P' j s at 
small Q^) and (iph describes the coupling of the Pomeron to a hadron, see Fig. ^ The 
normalization of the flux-factor jp is unfortunately ambiguous, since the Pomeron is 




t 



Figure 12. Schematic diagram of diffractive hadronic scattering via a partonic 
constituent of the Pomeron which assumes that the structure of the upper hadron 
can be factored off. an denotes the cross section of the 'hard' partonic subprocess. 
Unfortunately this partonic picture is not valid, since interactions with the remnants 
of the hadron cannot be neglected. 



not a real particle. This ambiguity poses serious problems in defining the correct 
normalization for the prediction of any hard cross section. According to section ^ we 
can define diffractive parton distributions for deep inelastic scattering. After dividing 
out the flux-factor one should arrive at a definition for the parton distribution of the 
Pomeron: 

jp(xp,t) 
fp{xp,t) 
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For real hadrons the momentum sum of all partons adds up to the total momentum of 
the hadron. This momentum sum rule provides an important normalization constraint 
on the determination of the parton distributions. Here there is no such constraint; ( p5| ) 
and ( ^6| ) only fix the relative quark versus gluon content in the Pomeron. Another 
problem arose when the predictons for diffractive processes were extrapolated from 
the CERN SPS-coUider energy {y/s = 630 GeV) to the Fermilab Tevatron proton- 
antiproton collider energy {^/s — 1800 GeV). The predicted cross sections overshot 
the measured values significantly. The origin of this disagreement has never been 
completely resolved. Unitarity corrections provide a possible explanation |9C| ]. For 
data from the HERA electron-proton collider, a consistent treatment of diffractive 
deep inelastic scattering and diffractive photoproduction of jets in terms of the parton 
distributions of the Pomeron seems possible [^if without renormalizing the flux- 
factor. However predictions for hard diffractive processes at the Tevatron, based on 
Pomeron distributions extracted from HERA data, fail badly unless one assumes a 
renormalization of the Pomeron flux- factor 

The problem is that although factorization has been proven for diffractive deep 
inelastic scattering, it is not valid in diffractive hadron scattering To discuss the 
issue of non-factorization in diffractive hadron scattering let us suppose that a quasi- 
real photon is a model for a hadron, and make use of our previous formalism. When 
we talk about factorization we have to specify three different cases: first, Regge-type 
factorization (|9^), second, coUinear factorization of the Pomeron and third coUinear 
factorization of the photon (hadron). We have already mentioned problems related 
to Regge-type factorization due to the ambiguity in the definition of the Pomeron 
flux-factor. In the genuine QCD-approach the Pomeron is solely viewed as a t channel 
vacuum exchange, and is not associated with real particles. The question whether or 
not a Pomeron trajectory exists, which passes through physical (glue-ball) states for 
positive values of t, remains open. 

A closer look at the diagrams presented in Fig. 5 shows the origin of the violation 
of coUinear factorization of the structure of the quasi-real photon, which we use as a 
model for the structure of the hadron. For example let us look at the triple Regge 
domain where we can assume a large transverse momentum in the final state while 
keeping the virtuality of the photon small and hadron-like. The left diagram of 
Fig . |l3| is one out of the complete set of diagrams which contribute to expression 
( p2| ) for diffractive qqg production. The particular feature is that only one quark line 
is involved in the scattering process, while the second has no interaction with a gluon. 
This second quark plays the role of the hadron-remnant. If we add up all related 
diagrams which leave the remnant quark unscattered, the gluon exchange amplitude 
simplifies from ( p3|) to 

f k, + k2 \ [dH f k2 + e , k2-e kA 

- U(fci + fe2)j. J ^ "^^(""'' ^ [jk^TW ^ (9^) 

where i,j = 1,2 are the vector components in the transverse plane. The first term, 
which we denote = (•■■)ii involves the quark propagator, and the second term 
bj = (• • ■)■ which is associated with the s channel gluon, is that shown in (|53| ) together 
with the t. — > — £ term. The dominant configuration is ki « — fe2, describing a quark 
and gluon jet in the final state, such that fcj ~ fc^ ^ (^1+^2)^ ~ Q^- The factorization 
in this case follows since 



(98) 
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Figure 13. The first diagram is an example of the subset of diffractive qqg 
production diagrams in which the t channel gluons couple to a single quark. A 
quasi-real photon is used to mimic an incoming hadron. It appears possible to factor 
off the 'photon' structure and to use the second diagram to describe the diffractive 
hadroproduction of two jets by means of the hard subprocess q + p ^ {qg) + p. The 
reason why this, in fact, is not possible is explained in the text. 



The first factor, originating from the = (• • in (|97|), leads to a logarithm in k\/Q'^ 
in the cross section. This factor would be divergent for — 0. It represents a 
collinear singularity which under usual circumstances is absorbed into the structure 
of the photon. 

It appears that exactly the same formalism is appropriate to describe diffractive 
hadron scattering with jets in the final state. It seems that we simply need to calculate 
the process g + p — *■ qg +p given by the second diagram of Fig. In the triple Regge 
limit, which in this case means 3> k^, one finds the same structure for the s 
channel gluon contribution as in the amplitude (j^j . This is of course expected due to 
the factorization of the photon structure mentioned earlier. By integrating over the 
azimuthal angle of £, after some algebra, we find: 

io 2n [{k,+£r + {k,-£r ki) - ^"^^ ^^^fcr ^ ' 

That is the low £ contribution cancels out such that the integration over i'^ has a 
lower hard cutoff given by fc|. The process q+p ^ qg +P thus seems to be calculable 
completely pertubatively. An identical conclusion has recently been made in Ref. |Q 
where the same process has been calculated without restricting it to the triple Regge 
limit. The problem is that there is a non-negligible rescattering of the partons with 
spectator partons due to 'soft' gluon exchange. Free quarks or free gluons do not 
exist. Quarks or gluons can only be part of the initial state for those processes 
where the factorization theorem is applicable. The true initial state is colourless 
(hadrons, photons, etc.) and single partons are always accompanied by a remnant 
which eventually takes part in the interaction. This is precisely what happens in 
the case of diffraction and is the reason why factorization fails for diffractive hadron 
scattering iQ. One of the t channel gluons interacts with the remnant, as depicted 
in Fig. and upsets the subtle cancellation that leads to the result (|9g|). Taking 
into account all possible interactions which involve the remnant brings us back to the 
complete expressions ( p2| ) and ( |5^ ) of Section ^ 

Insight into the problem may be obtained from a simple estimate of the cross 
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Figure 14. A contribution involving the scattering of the remnant vuhich not only 
ruins the possibility of using parton distributions of the Pomeron to describe hadronic 
diffractive scattering, but, even more, prevents us factoring off the structure of the 
quasi-real photon that is mimicking an incoming hadron. 



section. To obtain the diffractive hadronic cross section we have to integrate from 
Qq up to ki « where Qo denotes the typical hadronic scale. The dominant 
contribution comes from the lower end of the integration so we may put £^ <C ^2 in 
(p3|) and obtain 



M,, 



as T{xip,f') 



1 

^2 



(fc2)»(fc2)j 
^2 



If we further assume T '-^ \ /P then we find the diffractive cross section| 
da^ 1 



dtdAPdk^ 



M2 Q, 



2 1.4 

"-2 



(100) 



(101) 



On the other hand if we were to use (|99|) to obtain the cross section from the 
perturbative subprocess q + p qg + p then we would obtain 



da{q+p^ qg+p) 
dtdM'^dkl 



1 



M2 fc6 ^ 



(102) 



which is power-suppressed by Qo/^2 relative to the true answer of ( |10lD ||[ This 
suppression is due to the cancellations in the angular integration of ( |9g| ) and the 
consequent restriction of the integration over l"^. After connecting a gluon to the 
remnant, on the other hand, the integration over ranges down to the soft scale 
Qo boosting the cross section by the factor fc^/Qo- We conclude that we must use 
the full non- factorizing formula (^^. This example clearly demonstrates that the 
well known factorization of the cross section into universal parton distributions and 
a 'hard' partonic subprocess is violated for colourless two-gluon exchange. The same 
conclusion holds for multiple gluon exchange. These findings are of course unfortunate 
with regard to practical applications, as they do not provide any quantitative 
description of hard diffraction in hadron collisions since we need to know the hadronic 
analogue of the virtual photon wave function. It offers, however, a qualitative 

t The fej integration over the interval (fci + fc2)^ J; Qq h ring,s an extra factor of Qq which partially 



cancels the apparent 1/Qq behaviour which follows from (IOC) 

X A related discussion on diffractive bb-prodnction i n had ron collisions can be found in 



the results are not directly comparable to (101) and (102) since the soft interaction w 
in Fig. 3 of bsl is neglected. 
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ith the rc 
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explanation why the physics of diffraction in deep inelastic scattering is fundamentally 
different from diffraction in hadron scattering. The best way to proceed towards a 
quantitative description of diffractive hadronic processes is to model the wave function 
of the incoming hadrons. 

Issues related to non-factorization are the rapidity gap survival probability 
| p6| , |97[ | in double diffractive dissocation, i.e. both hadrons dissociate, and unitarity 
corrections. Multiple gluon exchanges lead to absorptive corrections which limit the 
growth of the cross section with increasing energy. Once there are a lot of t channel 
gluon exchanges the probability of radiation into the gap is, of course, high as well. 
These soft raditions also reduce the fraction of high-p( jets with a rapidity gap in 
between them 

16. Conclusions 

If we assume that QCD is indeed the theory of strong interaction then the Pomeron 
itself is to be constructed from QCD. The simplest model to produce a colourless 
exchange is two gluon exchange and we have spun our discussion of diffractive processes 
around this model. It provides a good starting point for further sophistications such as 
gluon ladders and multiple gluon exchange. Many features in diffractive deep inelastic 
scattering such as the /3-spectrum and the Q^-dependence directly follow from the 
forms of the wave functions of the photon. In lowest order in perturbation theory the 
photon dissociates into a quark-antiquark pair (qq). The corresponding wave function 
is mainly built from the photon coupling to the quarks and the energy denominator 
of the virtual qq state. At higher order the wave function formalism can be extended 
to also include a gluon (that is the qqg state). Gluon radiation becomes important 
and even the dominant contribution for large diffractive masses M. Longitudinally 
polarized photons play a special role as they dominate the large /3-regime despite 
the fact that these contributions are of higher twist nature. The remaining regime 
of medium f3 is filled with the contribution from transversely produced qq-pairs as 
illustrated in Fig. ^. We have also discussed the triple Regge limit where the mass M 
of the diffractively produced system is much bigger than Q. This regime is interesting 
as it allows us to interpret the qqg final state as a system of two quark dipoles (the 
gluon is decomposed in colour space into a quark-antiquark pair) scattering on the 
proton. In addition it is of use in the analysis of diffractive 3-jet data. 

We have also discussed the issue of nonperturbative effects in open diffractive deep 
inelastic scattering. Multiple gluon exchanges have to be considered as corrections. 
It has been argued that these corrections only slightly affect the /3-distribution as the 
photon wave functions do not change. Up to colour factors multiple gluon exchange 
can be recast into the unintegrated gluon distribution. 

The photon wave function formalism can be straightforwardly carried across 
to vector meson production. Light vector mesons show a fairly flat ratio of the 
longitudinal versus transverse cross section at large which is a consequence of 
7 decreasing with Q^, see (|85|). The large scale, together with the projection on an 
exclusive final state, allows perturbative evolution to take place inside the Pomeron, 
i.e. the Pomeron can be associated with perturbative gluon ladders (the so-called 
hard or QCD Pomeron). As a consequence the energy dependence is much steeper 
than expected from a soft Pomeron. A similar conclusion can be drawn for J/^ 
photoproduction and open charm production. In this case the hard scale is given by 
the charm mass. Precise data are becoming available at HERA for these processes. 
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So far the measurements show evidence of the QCD expectations. Another closely 
related process is the photoproduction of photons or vector mesons at large momentum 
transfer \t\. In this scenario the large scale is provided by t, but unlike deep inelastic 
scattering t stretches along the whole gluon ladder. Therefore a different evolution 
scheme following the BFKL equation is at work. 

We have demonstrated that diffractive deep inelastic scattering is consistent 
with collinear factorization and we have introduced diffractive parton distributions 
accordingly. On the other hand we have argued against general parton distributions 
of the Pomeron because factorization fails in diffractive hadron-hadron collisions. In 
particular we have shown the importance of the remnant of the dissociating hadron 
in estimating the diffractive cross section. 

In summary, the variety of diffractive processes in DIS, and the number of different 
kinematic distributions that are accessible to experiment, offer a sensitive probe of 
QCD. The optical theorem means that the diffractive amplitude is proportional to 
the inclusive cross section. The dependence on the square of a quantity magnifies the 
sensitivity of the diffractive probe. Examples are 7*p Vp which depends on the 
square of the (skewed) gluon distribution of the proton, and jp J/ipp at large \t\ 
which magnifies the BFKL- type enhancement effect. Although much theoretical work 
remains to be done (such as the computation of higher order corrections), it is already 
clear that the accumulation of precise data for DIS diffractive processes will greatly 
deepen our understanding of QCD. 

17. Acknowledgements 

We thank Markus Diehl, Misha Ryskin and Thomas Teubner for reading through the 
manuscript and for numerous helpful comments. One of us (MW) thanks the UK 
Particle Physics and Astronomy Research Council for support. 

References 

[1] A. Donnachie, P.V. Landshoflf, Nucl. Phys. B231 (1984) 189. 

[2] ZEUS CoUab., Phys. Lett. B315 (1993) 481; B332 (1994) 228; Z. Phys. C68 (1995) 569; Eur. 
Phys. J. CI (1998) 81; C6 (1999) 43; 

HI CoUab., Nucl. Phys. B429 (1994) 477; Phys. Lett. B348 (1995) 681; NucL Phys. B472 
(1996) 3; Z. Phys. C69 (1995) 27; C75 (1997) 607; C76 (1997) 613. 
[3] CDF CoUab., F. Abe et aL, Phys. Rev. D50 (1994) 5550. 

[4] CDF CoUab., F. Abe et al., Phys. Rev. Lett. 78 (1997) 2698; 79 (1997) 2636; 80 (1998) 1156; 

DO CoUab., S. Abachi et aL, Phys. Rev. Lett. 72 (1994) 2332; 76 (1996) 734. 
[5] P.D.B. CoUins, Regge theory and high energy physics, Cambridge University Press, 1977; 

A. B. Kaidalov, Phys. Rep. 50 (1979) 157. 
[6] LY. Pomeranchuk, Sov. Phys. JETP 7 (1958) 499. 
[7] A.H. Mueller, Phys. Rev. D2 (1970) 1963. 

[8] A. Donnachie, P.V. Landshoflf, Phys. Lett. 191B (1987) 309, Erratum-ibid. 198B (1987) 590. 

[9] M. Dichl, Eur. Phys. J. C6 (1999) 503. 
[10] J.D. Bjorkcn, Lecture Notes in Physics 56, Springer Verlag 1975, p. 93. 
[11] F.E. Low, Phys. Rev. D12 (1975) 163; 

S. Nussinov, Phys. Rev. Lett. 34 (1976) 1286. 
[12] J. Gunion, D.E. Soper, Phys. Rev. D15 (1977) 2617; 

E.M. Levin, M.G. Ryskin, Sov. .L Nucl. Phys. 34 (1981) 619. 
[13] A.H. Mueller, Eur. Phys. ,1. Al (1998) 19. 
[14] S.J. Brodsky, G. Lepage, Phys. Rev. D22 (1980) 2157. 
[15] N. Nikolaev, B.C. Zakharov, Z. Phys. C49 (1990) 607. 
[16] N. Nikolaev, B.G. Zakharov, Z. Phys. C53 (1992). 
[17] J.D. Bjorken, J. Kogut, D. Soper, Phys. Rev.D3 (1971) 1382. 



39 



[18] M.G. Ryskin, Sov. J. Nucl. Phys. 52 (1990) 529. 

[19] A.H. Mueller, Nucl. Phys. B335 (1990) 115. 

[20] D.Yu. Ivanov, M. Wiisthofr, Eur. Phys. J. C8 (1999) 107. 

[21] M. Wiisthofr, Phys. Rev. D56 (1997) 4311. 

[22] K. Golec-Biernat, M. Wiisthoff, preprint DTP-99-20, March 1999. 

[23] S. Catani, M. Ciafaloni, F. Hautmann, Phys. Lett. B242 (1990) 97; Nucl. Phys. B366 (1991) 

135; S. Catani, F. Hautmann, Nucl. Phys. B427 (1994) 475; 

J.C. CoUins, R.K. Ellis, Nucl. Phys. B360 (1991) 3. 
[24] J. Bartels, J. EUis, H. Kowalski, M. Wiisthoff, Eur. Phys. J. C7 (1999) 443. 
[25] P.V. Landshoff, O. Nachtmann, Z. Phys. C35 (1987) 405. 
[26] HI Collab., C. Adloff et al., Z. Phys. C76 (1997) 613. 
[27] M. Diehl, Z. Phys. C66 (1995) 181. 
[28] A.H. Mueller, B. Patel, Nucl. Phys. B425 (1994) 471. 
[29] N. Nikolaev, B. G. Zakharov, JETP 78 (1994) 598. 
[30] N. Nikolaev, B.G. Zakharov, Z. Phys. C64 (1994) 631. 
[31] A. Bialas, R. Peschanski, Phys. Lett. B378 (1996) 302. 
[32] A. Bialas, H. Navalet, R. Peschanski, Eur. Phys. J. C8 (1999) 643. 
[33] J.C. Collins, Phys. Rev. D57 (1998) 3051. 
[34] A. Berera, D.E. Soper, Phys. Rev. D53 (1996) 6162. 

[35] W. BuchmuUer, A. Hebecker, Nucl. Phys. B476 (1996) 203. 

[36] F. Hautmann, Z. Kunszt, D.E. Soper, Phys. Rev. Lett. 81 (1998) 3333; hep-ph/ 9906284 

[37] J. Bartels, C. Ewerz, H. Letter, M. Wiisthoff, Phys. Lett. B386 (1996) 389. 

[38] M. Diehl, Z. Phys. C76 (1997) 49 9: Eur. Phvs. J. C4 (1998) 497. 

[39] J. Bartels, H. Jung, M. Wiisthoff, hep-ph/9903265, preprint DTP/99/10. 

[40] J. Bartels, M. Wiisthoff, Z. Phys. <J66 (1995) 157. 

[41] E.A. Kuraev, L.N. Lipatov, V.S. Fadin, Phys. Lett. B60 (1975) 50; Sov. Phys. JETP 44 (1976) 
443; Sov. Phys. JETP 45 (1977) 199; 

Ya.Ya. Balitsky, L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 
[42] Yu. Dokshitzer, Sov. Phys. 46 (1977) 641; 

V.N.Gribov, L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438, 675; 

G. Altarelli, G. Parisi, Nucl. Phys. B126 (1977) 298. 
[43] J. Bartels, H. Lotter, M. Wiisthoff, Z. Phys. C68 (1995) 121. 
[44] C. Ewerz, PhD-thesis, preprint DESY-THESIS-1998-025. 
[45] W. BuchmuUer, T. Gehrmann, A. Hebecker, Nucl. Phys. B537 (1999) 477. 
[46] A. Hebecker, preprint HD-THEP-99-12, May 1999. 
[47] A.H. Mueller, private communication. 

[48] S. Munier, R. Peschanski, Ch. Royon, Nucl. Phys. B534 (1998) 297. 
[49] E. Witten, Nucl. Phys. B104 (1976) 445. 

[50] E.M. Levin, A.D. Martin, M.G. Ryskin, T. Teubner, Z. Phys. C74 (1997) 671. 
[51] H. Lotter, Phys. Lett. B406 (1997) 171. 
[52] M. Diehl, Eur. Phys. J. CI (1998) 293. 

[53] M. Genovese, N.N. Nikolaev, B.G. Zakharov, Phys. Lett. B378 (1996) 347. 

[54] E. Levin, M. Wiisthoff, Phys. Rev. D50 (1994) 4306. 

[55] J.C. Collins, L. Frankfurt, M. Strikman, Phys. Rev. D56 (1997) 2982. 

[56] P.A.M. Guichon and M. Vanderhaeghen, Prog, in Part, and Nucl. Phys. 41 (1998) 125. 

[57] X. Ji, J. Phys. G24 (1998) 1181. 

[58] A.V. Radyushkin, Phys. Rev. D56 (1997) 5524. 

[59] A.D. Martin, M.G. Ryskin, Phys. Rev. D57 (1998) 6692. 

[60] K.J. Golec-Biernat, A.D. Martin, Phys. Rev. D59 (1999) 014029. 

[61] A.G. Shuvaev, K.J. Golec-Biernat, A.D. Martin, M.G. Ryskin, Phys. Rev. D60 (1999) 014015. 
[62] J.R. Forshaw, M.G. Ryskin, Z. Phys. C68 (1995) 137. 

[63] J. Bartels, J.R. Forshaw, H. Lotter, M. Wiisthoff, Phys. Lett. B375 (1996) 301. 
[64] I.F. Ginzburg, D.Yu. Ivanov, Phys. Rev. D54 (1996) 5523. 

[65] L. Frankfurt, A. Freund, M. Strikman, Phys. Rev. D58 (1998) 114001; erratum, ibid D59 (1999) 
119901. 

[66] A. Donnachie, P.V. Landshoff, Phys. Lett. B185 (1987) 403; Nucl. Phys. B311 (1989) 509; 

Phys. Lett. B348 (1995) 213. 
[67] J.R. Cudell, Nucl. Phys. B336 (1990) 1; J.R. Cudell, I. Royen, Phys. Lett. B397 (1997) 317; 

G. Kerley G. Shaw, Phys. Rev. D56 (1997) 7291; 

D. Schildknecht, G.A. Schuler, B. Surrow, Phys. Lett. B449 (1999) 328. 
[68] E.L. Berger, D. Jones, Phys. Rev. D27 (1981) 1521. 



40 

[69] M.G. Ryskin, Z. Phys. C57 (1993) 89. 

[70] S.J. Brodsky, L. Frankfurt, J.F. Gunion, A.H. Mueller, and M. Strikman, Phys. Rev. D50 (1994) 
3134. 

[71] L. Frankfurt, W. Koepf, M. Strikman, Phys Rev. D54 (1996) 319; ibid. D57 (1998) 513. 

[72] M.G. Ryskin, R.G. Roberts, A.D. Martin, E.M. Levin, Z. Phys. C76 (1997) 231. 

[73] P. Hoodbhoy Phys. Rev. D56 (1997) 388. 

[74] A.D. Martin, M.G. Ryskin, T. Teubner, Phys. Rev. D55 (1997) 4329. 

[75] A.D. Martin, M.G. Ryskin, T. Teubner, Durham preprint DTP/99/98. 

[76] A.D. Martin, M.G. Ryskin, T. Teubner, Phys. Lett. B454 (1999) 339. 

[77] L. Frankfurt, M. McDermott, M. Strikman, JHEP (1999) 9902:002. 

[78] D.Yu. Ivanov, R. Kirschner, Phys. Rev. D58 (1998) 114026. 

[79] K. Schilling, G. Wolf, Nucl. Phys. B61 (1973) 381. 

[80] E.V. Kuraev, N.N. Nikolaev, E.G. Zakharov, JETP Lett. 68 (1998) 696. 

[81] L.N. Lipatov, Sov. Phys. JETP 63 (1986) 904. 

[82] V.S. Fadin, L.N. Lipatov, Phys. Lett. B429 (1998) 127. 

[83] M. Ciafaloni, G. Camici, Phys. Lett. B430 (1998) 349. 

[84] V.S. Fadin, R. Fiore, Phys. Lett. B440 (1998) 359. 

[85] I.F. Ginzburg, D.Yu. Ivanov, Phys. Rev. D54 (1996) 5523. 

[86] A.H. Mueller, W.K. Tang, Phys. Lett. B284 (1992) 123. 

[87] J. Bartels, J.R. Forshaw, H. Letter, L.N. Lipatov, M.G. Ryskin, M. Wiisthofl, Phys. Lett. B348 
(1995) 589. 

N.G. Evanson, J.R. Forshaw, \iep-ph/9902481 



[89] G. Ingelman, P.E. Schlein, Phys. Lett. 152B (1985) 256. 
[90] E. Gotsman, E.M. Levin, U. Maor, Phys. Rev. D49 (1994) 4321. 
[91] L. Alvero, J G Gnllins J Terr nn J.J. Whitmore, Phys. Rev. D59 (1999) 074022. 



[92] K. Goulianos, \iep-ph/9502S5f{ 
[93] F. Yuan, K. Chao, \iep-ph/9^nj.9,^\ 



[94] D.E. Soper, Proc. of DIS97, hep-ph/9707384 , 



[95] G. Alves, E.M. Levin, A. Santoro, Phys. Rev. D55 (1997) 2683. 

[96] J.D. Bjorken, Int. J. Mod. Phys. A7 (1992) 4189; Phys. Rev. D47 (1993) 101. 

[97] E. Gotsman, E. Levin, U. Maor, Phys. Lett. B438 (1998) 229. 



